Observability of Earth-skimming Ultra-high Energy Neutrinos 
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Neutrinos with energies above 10* GeV are expected from cosmic ray interactions with the mi- 
crowave background and are predicted in many speculative models. Such energetic neutrinos are 
difficult to detect, as they are shadowed by the Earth, but rarely interact in the atmosphere. Here 
we propose a novel detection strategy: Earth-skimming neutrinos convert to charged leptons that 
escape the Earth, and these leptons are detected in ground level fluorescence detectors. With the 
existing HiRes detector, neutrinos from some proposed sources are marginally detectable, and im- 
provements of two orders of magnitude are possible at the proposed Telescope Array. 
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Cosmic neutrinos with energies above 10® GeV, so 
far unobserved, have great potential as probes of as- 
trophysics and particle physics phenomena. They es- 
cape from dense regions of matter and point back to 
their sources, thereby providing a unique window into 
the most violent events in the universe. Once they reach 
the Earth, they interact with center-of-mass energies far 
beyond foreseeable man-made colliders and probe new 
physics at and beyond the weak scale. 

The sources of ultra-high energy neutrinos range from 
the well-established to the highly speculative The 
cosmic ray spectrum is well-measured up to the Greisen- 
Zatsepin-Kuz'min (GZK) cutoff [|| at 5 x IQi" GeV. 
Such cosmic rays necessarily interact with the 2.7°K cos- 
mic microwave background through pion photoproduc- 



tion p7 



, producing "Greisen neutrinos" when the 



pions decay In addition to this 'guaranteed' flux, 

far larger fluxes are predicted in models of active galactic 
nuclei (AGN) |^,^ and in proposed explanations of the 
observed cosmic rays with energies above the GZK cutoff. 
The latter include decays of topological defects (TDs) 
and Z- bursts ||] . Fluxes from photoproduction and some 
representative hypothesized sources are given in Fig. ^ 

The detection of ultra-high energy cosmic neutrinos is, 
however, extremely difficult, especially for those with en- 
ergies above 10® GeV. At these energies, the neutrino 
interaction length is below 2000 km water equivalent in 
rock, and so upward-going neutrinos are typically blocked 
by the Earth. This shadowing severely restricts rates in 
underground detectors such as AMANDA/IceCube 
which are bounded by detection volumes of at most 
1 km^. At the same time, the atmosphere is nearly 
transparent to these neutrinos. Even for quasi- horizontal 
neutrinos, which traverse an atmospheric depth of up to 
360 m water equivalent, fewer than 1 in 10^ convert to 
charged leptons. At the Pierre Auger Observatory, esti- 
mated detection rates are of order 0.1 to 1 event per year 
for Greisen neutrinos from the ground array 13 1, with 
similar rates for the Auger air fluorescence detectors . 

Here we explore an alternative method for detecting 
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FIG. 1. Differential fluxes of muon neutrinos (i/^ -I- i/^) 
from Greisen photoproduction (solid), active galactic nu- 
clei Q (long dashed), topological defects (short dashed), 
and Z-bursts (dotted). For maximal v^~Vr mixing, these 
fluxes are divided equally between ^ and r neutrinos when 
they reach the Earth. 



il2 



ultra-high energy neutrinos with > 10® GeV. While 
upward-going neutrinos are usually blocked by the Earth, 
those that skim the Earth, traveling at low angles along 
chords with lengths of order their interaction length, are 
not. Some of these neutrinos will convert to charged lep- 
tons. In particular, muon and tau leptons travel up to 
0(10 km) in the Earth at these energies, and so a signifi- 
cant number of them may exit the Earth and be detected 
by surface fluorescence detectors. A schematic picture of 
the events we are considering is given in Fig. ^. This 
method exploits both the Earth as a large-volume con- 
verter, and the atmosphere as a large-volume detector. 

Upward-going air showers have been discussed pre- 
viously [Tsl] , with an emphasis on differences between 
upward-going and conventional showers and the possibil- 
ity of space-based detection. The question of rates was 
not addressed. Very recently, the detection of showers 
from r decays in the Auger Observatory ground array 
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FIG. 2. A neutrino v enters the Earth with energy Ei, at 
nadir angle 9 and azimuthal angle (j). It then travels for dis- 
tance z before converting to a charged lepton £, which exits 
the Earth with energy Ei. 

has been considered The possibiUty of detecting 

moon-skimming neutrinos through radio signals has also 
recently generated interest [ p7||T^ ]. 

Given an isotropic neutrino flux the resulting dif- 
ferential flux of charged leptons exiting the Earth is 



d^e{Ei, cos 0,1 
dEp d cos 6 da 



where K is the probability that a neutrino entering the 
Earth with energy and nadir angle 9 produces a lep- 
ton that exits the Earth with energy Eg. Such an event 
requires that (a) the neutrino survives for some distance z 
in the Earth, (b) the neutrino then converts to a lepton, 
(c) the created lepton exits the Earth before decaying, 
and (d) the lepton's energy and position when produced 
are such that it leaves the Earth with energy Ee. 

The probability for a neutrino with energy Ei, and 
nadir angle 9 to survive for a distance z is 



Pa 



exp 



dz' 



(2) 



where L''cciEu,9, z) = [a!^ciE.)piri9, z))Na]-^ is the 
charged current interaction length, with a^Q^E^) the 
interaction cross section a{i>N £X) for a neutrino 
with energy Ei,, p{r) the Earth's density at distance r 
from its center, and Na = 6.022 x 10^3 g-i. The dis- 
tance r is given by r^{9,z) = i?^ + z^ ~ 2Rq)Z cos9, 
where « 6371 km is the radius of the Earth. For 
Ev ^ 10^ GeV, the charged current v and D cross sec- 
tions are virtually identical, and we may neglect multiple 
charged current interactions and neutrino energy degra- 
dation from neutral current processes. Also, at these 
energies, the optimal nadir angle for charged lepton pro- 
duction is 90° —9 « 1°. Leptons produced by Earth- 
skimming neutrinos travel essentially horizontally. 

The probability for neutrino conversion to a charged 
lepton in the interval [z, z+dz] is dz/ LQf~,{E^, 6, z). How- 
ever, given that detectable leptons travel nearly horizon- 
tally with path length of 0(10 km), this conversion must 
take place near the Earth's surface where the Earth's 



density is ps = 2.65 g/cm^. The conversion probability 
is then well-approximated by 



Ph 



dz 



where L 



CCs 



'cc 



{E,)psNa]- 



(3) 



We assume the 



lepton takes all of the neutrino energy. For ultra-high 
energy neutrinos, the mean inelasticity parameter is (1 — 
Ee/E^) « 0.2 ||l^. We therefore expect this assumption 
to make only a small difference. 

The survival probability Pc for a charged lepton losing 
energy as it moves through the Earth is described by the 
coupled differential equations 



dEe/dz = -{ae+ PeEi) p{r{9, z)) 
dPc/dz = -Pc/{cTiEe/mi) , 



(4) 
(5) 



where c is the speed of light, and mi and T£ are the 
lepton's rest mass and lifetime, respectively. Equa- 
tion (Q) parameterizes lepton energy loss through 
bremsstrahlung, pair production, and photonuclear in- 
teractions, under the assumption of uniform energy loss. 
For the energies of interest here, /J^- « 0.8 x 10^® cm^/g, 
« 6.0 X 10^6 cmVg |20[|2|], and the effects of a^,^ 
are negligible. At the Earth's surface, taus and muons 
lose a decade of energy in 11 km and 1.5 km, respec- 
tively. These differential equations are easily solved for a 
constant density ps, and the survival probability is 



Pc = exp 



mp 



CTlPtPs 



1 



(6) 



Muon lifetimes are long enough that Pc — 1, but this 
factor may play a significant role for taus. 

Finally, the lepton's energy and location when pro- 
duced must be consistent with an exit energy Ep. From 
Eq. (^, for constant density ps and negligible at, this 
condition is enforced with the delta function 



Combining Eqs. (1), (|), (|), and (0), the kernel is then 

K(E„9;E,)= I 
Jo 



r^Rffi cos 6 



PaPbPcPd 



(8) 



However, Eq. (^) may be further simplified, because the 
lepton's range in Earth is far less than the typical neu- 
trino interaction length. The kernel is therefore domi- 
nated by the contribution from z w 2i?0 cos 9, and we 
may replace z with 2i?0 cos 9 in Pa- The only remaining 
z-dependence is in Pd- Using J dzS{h{z)) = \dh/dz 
the z integration yields 



1 

h=0' 



K{E„9;E,) 



1 



mi 



X exp 



Jo 


cos e 


dz' 




(E^,»,^') 






1 






EpfiiPs 



(9) 



2 



In our calculations, we use the kernel of Eq. (g) with 
the Preliminary Earth density profile P4 I . Our cross sec- 
tion evaluation closely follows Refs. p^j2^ ; details will 
be presented elsewhere p^ ]. 

Muons and taus may be detected in fluorescence de- 
tectors either directly or indirectly through their decay 
products. We have evaluated rates for all of these possi- 
bilities p^ ; here we concentrate on the most promising 
signal from electromagnetic energy in r-decay showers. 
The recent discovery of near-maximal v^-v^ mixing | |2^ ] 
implies that, even at the high energies of interest here, 
Vt — ^ '■ ^ the Earth's surface. We assume also 
that tau decay initiates an electromagnetic shower with 
probability Bem = 80% and a typical energy, averaged 
over all r decay modes weighted by branching fraction, of 
-E^EM = \Et- At energies of order 10^° GeV, the typical 
shower length is ~ 10 km in the low atmosphere. 

We follow the analysis of Ref. to estimate the ef- 
fective aperture for r-decay induced showers. The signal 
from an electromagnetic shower must compete with the 
average noise from the night sky. By considering the 
signal to background ratio in individual photomultiplier 
tubes, the energy required for an electromagnetic shower 
to be detected was found to be 



10^ 



(10) 



where Rp is the shower's impact parameter in km, and 
\r « 18 km is the Rayleigh scattering length Ed 
is an energy characteristic of the detector. In particular, 
Ed oc s/MJiy^ cx ^/d/D^, where Ae = d/D is the angu- 
lar acceptance of each photomultiplier tube, and d and D 
are the diameters of the photomultiplier tubes and mir- 
ror aperture, respectively. For Fly's Eye, requiring a 4ct 
triggering threshold, the value Ed = 10^ GeV was ver- 
ified to reproduce the experimental data well |^^. For 
HiRes, D has been increased from 1.575 m to 2.0 m and 
d reduced from 14.4 cm to 3.5 cm we therefore take 
Ed « 3.2 X 10'' GeV. For each module of the proposed 
Telescope Array, and requiring a 4(t signal, Ed has been 
estimated to be roughly 4 x 10^ GeV Finally, the 

fluorescence detectors of the Auger Observatory will 
also be sensitive to Earth-skimming events; we expect 
their sensitivity to lie somewhere between that of HiRes 
and Telescope Array. 

Following Ref. |Q, we assume that showers are de- 
tected if and only if initiated within distance Rp of the 
detector. We also make use of the fact that, at these ener- 
gies, all r leptons exit the Earth horizontally. Apertures 
for Earth-skimming taus for each of the three detectors 
discussed above are given in Fig. ^. In each case, the 
aperture rises with energy until time dilation causes taus 
to decay too late to be detected. The HiRes aperture 
peaks at 2000 km^ sr near 3 x 10^° GeV. With increased 
sensitivity, however, the aperture peak rises and moves 
to lower energies, significantly enhancing detection rates. 
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FIG. 3. Effective aperture estimates for the detection of 
Earth-skimming r leptons through their decays to electromag- 
netic showers for Fly's Eye (dotted), HiRes (one site) (solid), 
and Telescope Array (one site) (dashed). 
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Given the kernel function K{E^,0; Er) and effective 
apertures (Ar2)off (-Er), the number of tau leptons de- 
tected is Nr ^ J dE^dErd cos ed(t) cose !^K{An)cgTD, 
where d^,^{E^)/dE^ is the differential flux originating 
from a given neutrino source, T is the time an experi- 
ment runs, and D is the duty cycle. To account for the 
requirement of clear moonless nights for fluorescence de- 
tection, we take D — 10%, corresponding to an observing 
period of 3 x 10^ s per year. 

Event rates for the four neutrino sources given in 
Fig. |l|, binned by tau energy, are summarized in Table |. 
[Note that these rates are suppressed relative to those 
presented in an earlier version of this paper.] We assume 
2 and 11 detectors for HiRes and Telescope Array, re- 
spectively; some reduction from overlapping fields of view 
may be expected. For HiRes, we find that neutrinos from 
AGN and TDs are marginally detectable. For Telescope 
Array, these rates are enhanced by more than two orders 
of magnitude — several Greisen neutrinos per year can be 
detected, and tens to hundreds of AGN and TD neutri- 
nos are possible. Note that the rates may be significantly 
enhanced by including multi-bang events, which we have 
neglected, and also if the r energy loss, dominated by un- 
certain photo-nuclear interactions, is less than our con- 
servative assumption. 

Hundreds or even tens of events will shed light on many 
aspects of ultra-high energy astrophysics. The energy 
spectrum of detected events varies from source to source, 
as evident in Table |. With many events, the source en- 
ergy spectrum may be determined by deconvolving the 
observed spectrum with the kernel function. Note also 
that these rates may be improved with detectors that 
cover the sky densely very near the horizon or by fil- 
ters optimized for nearly horizontal events. Placement of 
detectors in valleys, which effectively enhances the con- 
version volume, may also improve detection rates. 
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TABLE I. Expected number of j^,— induced electromagnetic 
showers detected by atmospheric fluorescence. Three years of 
running with duty cycle D = 10% is assumed. 



Detector 


Er (GeV) 


Greisen 


AGN 


TD 




Flv's Eve 


10** — lO'^ 


nnsQ 
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Total 


0.026 


0.66 


0.20 


0.026 


Telescope 


l(f ~ Iff 


1.4 


230 


41 


1.0 


Array 


l(f - ioi° 


3.2 


50 


20 


2.3 




10" - 10" 


0.62 


0.93 


2.2 


0.93 




10" - 10" 






0.045 


0.094 




10" - 10" 






0.00035 


0.0033 




Total 


5.2 


280 


63 


4.3 



Earth-skimming neutrinos also open up other possibil- 
ities for detection. Cerenkov radiation provides an al- 
ternative signal for showers initiated by r decay. Con- 
ventional air shower arrays, which deploy a large num- 
ber of modules over a horizontal area, are not opti- 
mally adapted to Earth-skimming events. It is interest- 
ing to contemplate 'vertical' arrays, say on the side of a 
mountain, that would intercept Earth-skimming showers 
originating from a very large surrounding area. Earth- 
skimming events may also be detected from space, as in 
the OWL/Airwatch proposal 
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